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Abstract. Non- linear gravitational clustering in a universe dominated by dark 
energy, modelled by a 'quintessence' scalar field, and cold dark matter with space-time 
varying mass is studied. Models of this type, where the variable mass is induced by 
dependence on the scalar field, as suggested by string theory or extra-dimensions, have 
been proposed as a viable solution of the coincidence problem. A general framework for 
the study of the non-linear phases of structure formation in scalar field cosmologies is 
provided, starting from a general relativistic treatment of the combined dark matter- 
dark energy system. As a first application, the mildly non-linear evolution of dark 
matter perturbations is obtained by a straightforward extension of the Zel'dovich 
approximation. We argue that the dark energy fluctuations may play an active role 
in cosmological structure formation if the scalar field effective potential develops a 
temporary spinodal instability during the evolution. 
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1. Introduction 

There is increasing agreement on the picture of the Universe on large scales: the 
analysis of the magnitude-redshift relation for type la Supernovae pQ, measurements of 
the anisotropies of the Cosmic Microwave Background (CMB) j2], and analyses of the 
galaxy distribution in large redshift catalogs 3\ all converge to show that the Universe 
is spatially flat and composed by roughly 30% of cold dark matter (DM), responsible 
for the structure we see via the small amount of baryons, and, for the remaining 70% 
by an unclustered form of dark energy (DE), responsible for its present-day accelerated 
expansion. In order to alleviate or avoid fine-tuning of initial conditions, an alternative 
solution to the usual cosmological constant has been proposed in the form of a dynamical 
vacuum energy component with negative equation of state, the so-called quintessence 

The most common candidate for quintessence is a self-interacting scalar field </>. 
Different choices for the effective potential have been considered in the literature. In 
order to avoid fine tuning on the initial conditions, particular attention has been devoted 
to potentials leading to scaling or tracker behavior, like those depending on the scalar 
field exponentially jH] or with a negative power-law 

In addition to self-interactions, the DE scalar field might in principle be coupled 
to any other particle or field present in Nature. However, given the extremely low 
mass typically expected for today, ~ H ~ 10~ 33 eV, its couplings to common 
particles, such as baryons or photons, are severely bounded by experimental tests on 
the equivalence principle and on the time variation of coupling constants. On the other 
hand, the coupling to non-baryonic DM is not so tightly constrained and might play a 
significant role in the cosmological evolution [Jj. In this paper we will take into account 
such a possibility by allowing the DM particles' mass to depend on the DE scalar field. 
The DM particles are then Varying Mass Particles (VAMPS), a possibility that was 
discussed in the past in connection to the age problem of the Universe jS]. 

The VAMPS scenario was considered more recently in the quintessence context in 
jSj, where it was shown that it may lead to a solution of the so called 'cosmic coincidence' 
problem. Indeed, assuming an exponential form both for the coupling and for the 
effective potential, the late time cosmology may be such that the universe accelerates 
and the energy densities of DE and DM scale at the same rate. As expected, in this 
scenario the predictions for the DM abundance are drastically altered with respect to 
the standard freeze-out scenario jH] • A coupling of this type may arise for radii moduli in 
braneworld scenarios ^Oj- Couplings between DE and DM may also emerge in strongly 
coupled string theory [TT] . 

It is generally believed that perturbations of the quintessence scalar field play a 
negligible dynamical role in cosmic structure formation, because of the extreme smallness 
of the quintessence mass (m^ ~ H) which implies that its spatial fluctuations only 
appear on very large scales and are bound to be linear. Nevertheless, a number of 
models have been proposed where DE perturbations might in principle grow non-linear 
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and play an active role down to galactic scales. These include coupled-quintessence 
[T2~] , extended quintessence [T3J E] , inhomogeneous Chaplygin gas [T5| , k-essence [TH] . 
A key question in the context of scalar field cosmologies is whether fluctuations of the 
quintessence field can influence the galaxy formation process in a sensible way, affecting 
the dynamics of structure formation during the non-linear phases e.g. in the cores of 
dark matter halos. This problem has been recently addressed by Wetterich ^7], but 
the issue of whether quintessence fluctuations can grow non-linear on small scales and 
actively influence the non-linear clustering of the CDM component is still completely 
open. 

The VAMPS scenario considered here provides an useful bench mark to study 
situations where non-linearities in the scalar field may be induced either by the field's 
self interaction or by its couplings to other cosmic components. 

The plan of the paper is as follows. In Section 2 we give a full derivation of the 
equations which govern the non-linear evolution of DM and DE perturbations on sub- 
horizon scales. Then, in order to gain some intuitive insight on the content of these 
equations, we specialize to the VAMP model discussed in . A simple solution of the 
equations in the mildly non-linear regime for this model is obtained in Section 3, in 
terms of first-order lagrangian perturbation theory. In Section 4 we discuss how the 
simple model considered here could be modified to give rise to non-negligible scalar field 
fluctuations on scales relevant for galaxy or galaxy-cluster formation. 

2. General formalism 

As a starting point, we will expand the metric tensor g^ u (greek indices run from to 

3, while latin ones from 1 to 3) to linear order around a flat Robertson- Walker (RW) 
metric, while keeping the full - i.e. non-linear - energy momentum tensor on the RHS of 
Einstein equations. In the conformal Newtonian gauge, vector and tensor perturbations 
of the metric decouple from the scalar ones as long as they are treated linearly; by taking 
the spatial covariant divergence of the — i Einstein equations and the trace of the % — j 
equations, we get rid of vector and tensor modes and single out scalar ones. Thus, it is 
not restrictive to consider the metric [TH] 



where the scalar perturbations $ and \I/ are related to the gravitational potential, and 
are of order |u| 2 <C 1 (u being the typical velocity of DM particles) even in the presence 
of highly non-linear DM overdensities. Here a(rj) is the background scale-factor and r\ 
is the conformal time. 

In this approximation, Einstein's equations read 



where dR^ is the linear perturbation to the Ricci tensor, while =T^ — \8^ V T is the 
fully non-linear source and Stf is the background one. 



ds 2 = a 2 (r]) (1 + 2$)dr/ 2 - (1 - 2V)6 ij dx i dx j 



(1) 



(2) 
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The Friedmann equations for the background RW metric read 

H 2 = (p m + P</>) (3) 

n = — a (p m + p$ + 3^0) , (4) 

where H = a/a, and the dot denotes differentiation with respect to 77. 

We will take into account the contribution to the energy-momentum tensor coming 
from the dark matter and the scalar field. The former is made up of a discrete set 
of particles with field-dependent mass 771(0(77, x)) and coordinates x a (?7) (a = 1,2, . . .), 
described by the energy-momentum tensor 

T» ~ a~ 2 p m (0) u»u v , (5) 

where 

p m (0) = m(0) a' 3 5(3) ( x - x fl ) , 

a 

= dx^/dr] (x° = 77), and we have neglected 0(|u| 2 ) corrections. 
The scalar field energy-momentum tensor reads 

( M 2 \ 
T% = M 2 ^ .„ - St l-f<f> ;a ;a - V (</>)) , (6) 

where is dimensionless because of the factorization of M 2 = (87rG) _1 . After replacing 
the background quantities g^ v , 0, and in © and © one can read out the background 
energy density and pressure for our two components, p m , p^, and p^. 

The Einstein equations provide a redundant set of equations for the evolution of 
the linearized metric perturbations; we choose to use the — equation and the spatial 
covariant derivative of the — z's. They read, respectively, 



a 2 

V 2 $ + 3ft($ + *) + 3^ = — - 

2M p 



2M 2 ^2 

Pm {<t>) - Pm + —f{ct> 2 - ) - 2(y(0) - V{4>)) 

a z 



HV 2 <S> + V 2 * 



a 2 
2M 2 



M 2 2M 2 
di (pmttW) + ^(^09> + 0V 2 0) + —5.^*9* 



, (7) 



(8) 



a* a' 

where primes denote differentiation with respect to and, from now on, latin indices 
are raised with a Kronecker delta. 

The evolution of the scalar field 0(77, x) is given by the Klein-Gordon equation, 
which, for |$|, \u l \ <C 1, reads 

+ 2^0 - V 2 + |^'(0) = V($ - *) • V0 , (9) 

p 

where the effective potential £7(0) is the sum of the scalar field potential and the field- 
dependent dark matter energy density, 

rj(0)=y(0) + Pm (0) (10) 

Finally, we need the equation of motion for the dark matter particles, 

d 2 x fl ^ ^ dx a dx@ m ' Q^ff) 
ds 2 CT/3 ds ds m 
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where the last term comes from the field dependence of the particle mass and makes the 
trajectories deviate from geodesies. The spatial coordinates of the a-th particle satisfy 

77?' 

= -<9'$ &<f>, (11) 

771 

where H, and (ft are evaluated at (77, x a (?7)). 

Equations 0, ©, © and (|TT| form a close system in the variables $, <fi and 

As we are interested in sub-horizon scales, there are further approximations which 
can be made. First of all, in Eq. (JZj), only the first term on the LHS will be kept, as the 
others are C(7i 2 $ , TC 2 ^). Furthermore, the RHS of Eq. © is suppressed by at least 
an factor with respect to the Laplacian on the LHS, and can then be dropped. 

As a result, the \1/ perturbation is not needed anymore, $ can be identified with the 
gravitational potential, and we can concentrate on the reduced system of Eqs. ©, 0, 
and (fTTjl in their approximated forms. 

To illustrate the role of the scalar field we split it into background, 'non-relativistic' 
and 'relativistic' contributions, 

0(7/, x) = 0(77) + 0Ar/?(r/, x) + (p R (r), x) . (12) 

0(77) is the solution of the background Klein-Gordon equation, 

| + 2ft0 + a 2 M-V(0) = 0, (13) 

while the non-relativistic part 4>nr{Vi x ) is defined as the solution of the Poisson-like 
equation 

V 2 <j) NR = a 2 M p 2 AU\4> + <p NR ) , (14) 

where AU(4>) = U(<p) — V(4>) — p m . Thus, 4>nr accounts for the response of the scalar 
field to the localized DM distribution, as induced by the field dependence of the DM 
particle mass. It behaves similarly to the gravitational potential $ and, as we will see, 
in specific cases it may turn out to be just proportional to the latter. 
The remaining 'relativistic' part solves the equation 

'<P R +2H^ R -V 2 <p R +'<p NR +2H^ NR +a 2 M; 2 [AU\^+<p NR +<p R ^ =0.(15) 

The effective mass for <fi R , as read out from Eq. (fTB]) is given by AU"((j) + 4>nr), 
which, for the scales of interest for the formation of cosmic structures is much smaller 
than the corresponding k 2 /a 2 . Moreover, since the time-scale for <f> and (f>NR variation is 
given by H, the 4>nr and (pNR terms in (fTK|) provide 0{7i 2 ) sources. As a consequence, 
the (j) R component behaves essentially as the superposition of approximately massless 
plane waves on the scales of interest, with 0(a 2 p m /k 2 M 2 ) -C 1 amplitudes, and can 
then be neglected. 

The relevant dynamics can then be described in terms of the equations 



m 
m 
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V 2 $ 



a 



2 



2M p 2 



Pm(0 + <Pn R ) -Pm-2 (V(4> + <j) NR ) - V(<f>))] , (16) 



V 2 <j> NR = -^MJ\4> + <i> NR ), (17) 



x\ + 



p 

771 — 

n + —{4> + 4>nr) 

m 



TH — 

-d^-—d\4> + 4> NR ), (18) 

m 



and the background ones, Eqs. (jSJ), ©, and 

The equations above extend the standard Newtonian approximation from the pure 
DM case to the coupled system of DM and quintessence DE. All non-background DE 
effects are given by <Pnr, which solves the Poisson-like eq. (JT7J). 

In order to get a flavour of the dynamical content of the system, we specialize to a 
simple model which was proposed as a solution of the coincidence problem (HIE]]. The 
scalar potential and the DM particle mass are both exponentially dependent on (f>, 

V(<j>)=Vexp(J3<f>), m = mexp(-A0), (19) 

with (3, A > 0. The model has a solution such that 

— — 3 „ 3 + A(A + ,3) 

logo, n+ = i-n m = -, 2 ^ , (20) 



a + /3 b ' 9 m (\ + py 



which is an attractor in field space if (3 > (—A + \/X 2 + 12)/2. From the moment the 
attractor is reached, the energy densities in DM and in DE evolve at a constant ratio 
depending only on A and j3, thus solving the cosmic coincidence problem. 
Indeed, the solution (|2T)|) implies 

p m = m((f))n m ~ p ~ a- 3(1+w) , with w = - j^j , (21) 

where n m is the DM background number density. The negative w leads, if w < —1/3, to 
accelerated expansion of the universe. The <fi dependence of the DM mass modifies the 
usual scaling a -3 of non relativistic matter. Since the mass increases with the expansion, 
the effect is analogous to that of a fluid with negative equation of state, though DM is, 
as in the standard scenario, a pressureless fluid made up by non-interacting particles. 

The impact on the CMB of this model is similar to the one considered in jT^]. It 
was shown there that in order to get results in agreement with the CMB anisotropy 
spectrum a time-dependent A should be invoked, such that the early time dependence 
is the standard CDM one, and the attractor behavior of Eq. ()20)1 starts only recently. 
Since we use the model only for illustrative purposes, we will stick to a constant A in 
the following. 

The exponential dependences on in A£7 and Eq. (fTTjl imply that also 4>nr = 
0(a 2 p m /k 2 Mp) <^ 1. Then, as dark matter fluctuations approach the non linear regime, 
they start dominating over those of the effective potential, i.e. 

AU ~ m(0) ( a~ 3 £ 5 (3) (x - x a ) - n m ) . (22) 
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In this approximation, the RHS of the Poisson's equation ()16)) and that of Eq. (J 17)) are 
proportional to each other, so that (p NR = — 2A$ = 0(|u| 2 ). 

Using the properties of the attractor solution, Eq. (ffilj) . we get the further reduced 
system, 

V 2 $ = ^5 ^ («~ 3 E ^ (3) (x - x a ) - n m ) , (23) 
xl + (l-3w)Hx i a = -(l + 2\ 2 )d l <S> (24) 

The coupling between DM and DE modifies the equations from the standard ones in two 
ways. It changes the background evolution through the w-terms, and it gives an extra 
contribution to the force between two DM particles, due to the exchange of a light <fi 
quantum. The last effect is responsible for the 2A 2 correction to the gravitational force 
in the LHS of Eq. 



3. Mildly non-linear evolution 

The set of equations (|2HJ), (j2H) provides the required extension of the standard 
Newtonian approximation (e.g. [20J), describing the non-linear evolution of matter 
perturbations, to the case where a quintessence scalar field component is present. As a 
preliminary application, we can study the effect of the scalar field on the DM evolution 
in the midly non-linear regime, using first-order Lagrangian perturbation theory. 

Standing to standard notation [21], one can look for a time variable r = t(t]) 
by which the equation of motion derived from (jl8|) does not contain the Hubble drag 
term. We can then write an equation for the displacement from the initial (Lagrangian) 
position q to the final (Eulerian) position x = q + S(r, q) subjected to a force per 
unit mass proportional to <9*$. After taking the divergence of the resulting equation of 
motion, one can use the Poisson Eq. (j2HJ) to obtain an equation for the deformation 
tensor Si j(r,q) = dS l /dqi and then reconstruct the density field p(r, x) = p(q)J _1 , 
where J = det(5jj + Sij) is the determinant of the Jacobian dx % /dq J . 

On the attractor (j2U|) . for t = a n and n = — 1(1 + 9w), the DM equation of motion 
is simplified to 

d 2 S (1 + 2A 2 1 ^ 

and the Poisson equation ()2^|) becomes 




V xV = ( ^ (1 ^ 2 A2) ) Ux(r, q)) = b(r)5 , (26) 

where 5 = 5 m (x(r, q),r) = [p m (r, x) — p m (T)]/Pm(T) is the density contrast. 

Expanding the resulting equation to first order in the displacement vector S and 
making the ansatz S(r, q) = /(r)S(q) we obtain 

^/(r)-6(r)/(r)=0, (27) 
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whose solutions read 



/(r) = 



P 



2 ± 2\ 



1 + 



24fi m (l + 2A 2 
(1 -9w) 2 



(2? 



We also get V q • S(q) = — <5 , where £ is the VAMPS density fluctuation linearly 
extrapolated to the present time. 
The matter density field follows: 



Pn 



2T,X 



nL [i + A,,(q)(i + ^) 

where Aj(q) are the eigenvalues of dSj/dq % and 

9w) 



-m+ 



m± 



\ 



1 + 



240J1 + 2A 2 ) 



(29) 



(30) 



;i -9w) 2 

(in agreement with ^2]): the effect of the scalar field and a non- vanishing coupling 
constant A is to anticipate (if m + > 1) or delay (if m + < 1) the formation of pancakes 
with respect to scenarios with only DM with constant mass (m + = 1), depending on 
the value of A and (3. Because the DM-DE coupling in our specific model affects the 
motion of DM particles only through a 2A 2 correction of the gravitational force, the 
eingenvectors of the deformation tensor dSj/dq 1 turn out to be exactly aligned with 
those of the gravitational force itself. As a consequence, no other scalar field effects 
appear on small scales. 



4. Conclusions 



In this paper we have set the stage for future studies of non-linear gravitational clustering 
in scalar field cosmologies, extending the Newtonian approximation to the case of a 
coupled DM-DE system. The non-linearity in the scalar field may emerge either as a 
result of its self interactions or by its coupling to DM. It obeys a Poisson-like equation, 
Eq. (j!7|) . quite similar to that for the gravitational potential <£>. 

In order to get some analytic insight, we have specialized to the VAMP scenario 
discussed in [9/ , in which the DM particle's mass depends exponentially on the DE field. 
In this particular model, we found that the scalar field perturbations remain linear inside 
the horizon, although the coupling plays a non-trivial dynamical role. 

The mildly non-linear regime has been tested here, using first-order Lagrangian 
perturbation theory. The DM-DE coupling weakly affects the gravitational field, acting 
on the onset of pancake formation without any further consequence on the spatial 
distribution of the DM density field on the scales of interest. 

The constant coupling A considered here is, however, disfavored by comparison 
with the CMB anisotropy spectrum. A time- or field-dependent A - smaller in the past 
and approximately constant in the recent epoch - would take care of this problem, as 
discussed in jTH]- in this case, our results would not change qualitatively. Namely, the 
scalar field fluctuations would still remain well inside the linear regime and the onset of 



Non-linear gravitational clustering in scalar field cosmologies 



9 



structure formation would be changed with respect to ordinary quintessence, though of 
a lesser extent than for a constant A. 

Scalar field non-linearities are, however, not precluded in the VAMP scenario. For 
instance, modifying the field dependence of the DM mass in Eq. (|19p. the full potential 
U(<f) may temporarily develop a spinodal instability (negative second derivative). 
During this time, perturbations at scales k/a < \J—U" grow exponentially and easily 
exit the linear regime. When the spinodal epoch ends, the system can be described 
by our equations (fTHjl . (fTTjl . (fTHJ) . assigning the developed non-linearities to the initial 
conditions for the non-relativistic field fluctuation <Pnr- 

The approach developed in this paper can be straightforwardly applied to the 
numerical study of models where the quintessence field fluctuates non-linearly, as those 
proposed in refs. [121 CHI d EES EEE1 E! • 
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